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1. INTRODUCTION 

Photovoltaic systems are used widely in employment, either solely or in confluence with other 
electrical power sources. Applications powered by photovoltaic (PV) systems include communication 
equipment, somehow powered telecommunications installations, remote monitoring, lighting, water 
pumping, and battery charging. The modeling of PV arrays and modules plays an essential role in the 
performance and productivity of PV systems. The extraction of solar cell or PV module parameters is critical 
for determining the performance of various jobs and carrying out their design and quality control. 

Several techniques and approaches have been offered to obtain the five physical parameters using 
the one diode model. These approaches are grouped into three classes: analytical, numerical, and 
evolutionary approaches [1]-[23]. Villalva et al. [5], which is among the most cited diving factories in this 
area, proposed a simple algorithm to extract the five parameters. This algorithm uses the series resistance as a 
duplication parameter to minimize the error between the calculated and measured peak power values. The 
disadvantage of this method is that it’s accurate near the maximum power point but inaccurate in other 
regions and uses a fixed value of ideality factor equal to 1.3. Cubas et al. [11] has proposed a coherent 
approach using four analytical expressions with some approximation to extract the four parameter values and 
an ideality factor equal to 1.3. Ma et al. [17] have presented an extraction algorithm based on bio-inspired. 
Furthermore, a study of two technics using an iterative algorithm and the Lambert function has been 
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suggested in work [18] to determine the five parameters of the PV modules under varying environmental 
conditions. 

New, Achouby et al. [21] provided an exact numerical approach to find the five physical 
parameters. This method is based on varying the ideality parameter and resolving a nonlinear system of four 
equations. This method requires appropriate initial guess for four physical parameters. Zaimi et al. [22] 
presented a method based on a combination of well-founded and numerical approaches. This technique needs 
two coherent initial values of series resistance and ideality factor. Stornelli et al. [23] suggested a new, 
five-parameter method. This simplified method allows one to find the optimal value of the ideality parameter 
and the shunt resistance. 

In this investigation, we present another technique based on a combination of iterative and analytical 
approaches to extract physical parameters of a single diode model of a PV generator. In the first step, we use 
a well-founded expression to calculate the series resistance Rs. In the spare step, we derive three analytic 
equations giving the parallel resistor Ryn, the saturation current Is, and the photo-current Ipn as functions of 
series resistance R, and ideality factor n using the linear least squares method. In the final step, we use the 
ideality factor as a parameter iteration to minimize the normalized root mean square error (NRMSE). This 
procedure has the advantage of using only one foremost value of an ideality parameter and allows us to 
reduce the number of unknowns to one. 


2. PROPOSED METHOD 
2.1. PV module modelling 

A PV module consists of several solar cells in series or parallel that transform solar irradiation into 
electrical current. To describe the PV module behavior, we use in our work the one diode model with five 
parameters; this model is the most used in the literature. The electrical equivalent circuit of the model is 
depicted in Figure 1. 


Figure 1. Equivalent circuit for one diode PV cell model 


The mathematical behavior of the PV panel connecting the output current and voltage and the five 
parameters is given by (1). 


I= Ip- l (exp ("=") - 1) - (1) 


n.NsVep Rsh 


Where Ipn is the photoelectric current, Is is the current of saturation on the diode, Rsn is the parallel resistance, 
n is the ideality factor, Rs is the series resistance, Ns is the number of solar cell connected in series 


(N;=1 for the solar cell), Vin is the thermal voltage given by: Vi, = = Where q is the electronic charge, kg 


is the Boltzman’s constant and T is the temperature in Kelvin. The equation (1) is implicit and requires 
numerical resolution. The explicit expression of this equation using the Lambert function is given by (2): 


I=- 


Is+tlyh G n.Ns V Iç.R V+Rs.Ustl yn 
2P sh__y 4 stth LambertW ( ss exp ( si(lstlpi) )) (2) 
1+Rs-Gsh 1+Rs-Gsh Rs N.Ns Vth (1+Rs-Gsh) Nn.Ns-Vth(1+Rs-Gsh) 


where Gsn is the parallel conductance given by 1/Rsn. 
2.2. Extraction method 
2.2.1. Analytical expression of Rs 
The evaluation of the (1) at the three remarkable points of I-V curve (V=Voc, I=0), (V=0, I=Isc) and 


(V=V mp, I=Imp) points, respectively, provide the equations: 
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0 = Ip = Is (exp (rv) ~ 1) — Fe © 
ke-m- em (i) ~1) 8 a 
Imp = lyn — 1s (exp e ae sto (5) 


According to [11] and using the three previous equations, R, can be expressed analytically: 


aVmp+B (Voc—Vmp) 


Re ep lnp (6) 
a and ß are given by: 

a = (Vap + (Imp — fee) Rena) In (eae s0) | = Vap — Rerolmp o) 
Rsho is the slope at the short circuit point: 

Rao =~ Flay. (8) 


2.2.2. Calculation of Is, Iph and Gsh 
In order to calculate the three parameters Iph, Is and Rsn as function of R, and n, we use a linear least 
square method to minimize the follow error expression: 


F= Er bea (Vex) = ba) = 


Dita (lpn — I; (exp (eatis) ~ 1) = Gsn(Viex + Rs. Lex) ~ a) í a 


NVth 


Tiex and Iica are the measured and the theoretical current. Vi ex is the experimental voltage of the PV module 
and N is the number of experimental points. The minimization of the F function requires the system 
resolution (10): 


ôli, ca 
m 2 Li n Clica — liex) = 0 
ôli ca 
Ast 2 Disi alg Clica = lex) =0 (10) 


Dlica 
aGe, a = =25 1 IGsn Clica = lex) =0 
After some mathematical operations, the previous system of equations become [19]: 
Ip, N =l; Di =1 EXP; 5 Gsh EL 1 Vshi = Ðtli lex 


—Ipn Die EXP; + 1, Diy EXP;? + Gor Dia Vsn EXP; = — EiL lex EXP; (11) 
—Ipn Dia Vou, + 1s Diha Vou, EXP: + Gop Dik Vong = — Lia liexVshi 


where EXP; and Vni are given by (12): 


EXP; = exp (“otBstiex) _ 1; Veni = Vi ex F Rs. li ex (12) 


n.NsVth 


2.2.3. Algorithm of ideality factor extraction 

In this algorithm, we use ideality factor n as variation parameter which varies between 1 and 2 for 
silicon material and solve the linear system of (10) using MATLAB to find Iph, Is and Rsn values. Then, we 
select n that minimize the NRMSE value: 
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ua 

SIN Ciim TF 

NRMSE = [aztal f (13) 
W 2i=1lim 


In Figure 2, we have detailed the extraction strategy of the five parameters using our proposed technique. The 


main steps are given by the following flowchart. 


- Load experimental curve I (V) 

- Input values of: V mps Imp Voc and Ies 

- Input value of R,,o , equation (8) 

- Get value of R, , equation (6) 

- Define increment value An and initialize n = 1 


- Solve the system of equations (11) 

- Get values of I,, I, and Rop 

- Calculate theoretical values of current I, equation (2) 
- Calculate NRMSE value 


n=n+An 


- Find the optimal value of n that minimize the NRMSE 
- Get final values of I, I, and Ry, 


= 


Figure 2. Proposed flowchart method 


3. RESULTS AND DISCUSSION 

We confirm the suggested approach through two experimental test cases: the RTC solar cell at 33 °C 
and 1000 W/m? [1] and the multicrystalline panel KC200GT operating at standard test conditions (STC) 
(T=25 °C and G=1000 W/m2), the experimental current-voltage curve are extracted from manufacturer’s 
datasheet [24]. In order to evaluate the accuracy of our method, we use the NRMSE using (13) and the 
absolute error given by (14): 


AE = |li exp — lical (14) 


In Table 1, we summarized the remarkable points values and Rsno for the solar cell and the PV panel. 


Table 1. Remarkable point’s values and Rsno 


RTC solar cell KC200GT PV panel 
Imp (A) 0.67 7.61 
Vinp (V) 0.46 26.3 
Ix (A) 0.76 8.21 
Vo (V) 0.57 32.9 
Ro (Q) 58.656 597 


Tables 2 and 3 present the obtained values of the five parameters using the our technic for the solar 
cell and PV module respectively. For comprehensive comparison, the proposed algorithm is evaluated 
against previously techniques. According to this two tables, our method has the lowest NRMSE value. In 
Figures 3(a) and 3(b), we plot the NRMSE values versus ideality factor n, for the two experimental cases. 
According to these figures, the optimal value of n which minimizes the NRMSE is equal to 1.55 for solar cell 
and equal to 1.1 for PV module. 
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Figures 4(a) and 4(b) present the measured and the simulated current-voltage I (V) curves obtained 
by using the Lambert-function model, for the two cases studies, solar cell and KC200GT. As it can be seen 
for both cases, the simulated characteristics are in excellent accordance with the experimental data. The 
benefits of the proposed method are that need only one initial guess. 


Table 2. Solar cell parameters at 33 °C and 1000 W/m? 
Parameters Method [3] Method [6] Method [25] Proposed method 


RQ) 0.037 0.036 0.0355 0.033 
n 1.45 1.48 1.4905 1.55 
L(A) 0.7611 0.7607 0.7611 0.7606 
IA) 2.422.107 3.267.107 3.514.107 6.271.107 
R,(Q) 42 60.24 45.0472 66.36 
NRMSE 0.078 0.068 0.0072 0.0026 


Table 3. PV module parameters at STC conditions 
Parameters Method [11] Method [5] Method [19] Method [23] __ Proposed method 


Rs(Q) 0.23 0.221 0.233 0.2185 0.23 
n 13 13 1.0758 1.1 i 
Te) 8.213 8.214 8.211 8.196 8.211 
I(A) 9.76.10 9.825.108 2.12.10 3.27.10? 3.143.109 
Ry (Q) 597.38 415.405 132.88 164.2 154.64 
NRMSE 0.0571 0.0564 0.0086 0.0087 0.0069 
0.07 F : : i z 
0.06 
0.05|- 
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Figure 4. Measured and simulated characteristics I (V) for (a) solar cell and (b) KC200GT 
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In order to evaluate the impact of temperature cell T and irradiance level G on the five parameters 
and the I-V characteristic, we use the following expressions [19]. 


n(T,G) = Nsrc (15) 
Rs (T, G) = Rs strc (16) 
Rs (T,G) = Rsrsrc (17) 
Ipp(T, G) = =. (Ipnsre + Ki (T — Tsrc)) (18) 


fl 4(T,G)—G nVoc(T,G) 
I(T, G) = a ( TTO )- 


Ns-n.V aT) 


(19) 


Where, Gstc and Tstc are the overall radiation ant temperature under STC conditions and Voc is given by the 
following analytical model [21]: 


G 
Voc(T, G) = Voesre + Ky. (T — Tsc) +0-Ng.Ven(T). In (=) (20) 


Where, Ky is the temperature coefficient of Voc and Voc, stc is the open-circuit voltage in STC conditions. 

Figure 5(a) displays the measured and the simulated I-V characteristic for the KC200GT panel. As 
can be clearly seen, the simulated curve of the aforementioned PV module is in very close accord with the 
measured data at T=25 °C and at various irradiation levels. In Figure 5(b), we display variations of absolute 
error curve for a temparture 25 °C and different irradiation levels. In the voltage variation interval [0, Vimpl, 
the theoretical and experimental I=f(V) curves are in good agreement and the absolute current error tends to 
cancel. When, at a voltage variation n the interval [Vmp, Voc], the current decreases rapidly and the AE 
increases slightly but still remains below 0.15 for the irradiadiance 1000 W/m? and 0.3 for the iradainces 
800 W/m2, 600 W/m?, 400 W/m? and 200 W/m?. The current has an exponential characteristic in the zone 
[Vimp, Voc] at an almost constant voltage. 


O Experimental data T=25°C 
__ Simulated data 


G=1000w/m? 0.4 
se 2 0-2-6 o—S 


T=25°C 


o—o—6- 


6F 
G=600w/m? 
o-6—9—9—_9—9—_9—_90-9-9 


4+ 


I(A) 


G=400w/m? 
o—O—o- B-S-9—9—9-9— 0. 


2L G=200w/m? 
oo-e-86—-8—-8—-8_8 99.90. Sh 
40 
0 i 
0 5 10 15 20 25 30 35 
v(V) 
(a) 


Figure 5. Accuracy of the proposed method (a) measured and simulated I(V) characteristics and 
(b) absolute error AE(V), for KC200GT panel under varying irradiance 


In Figure 6(a), we display the measured and the simulated I-V characteristic, for the KC200GT 
panel for fixed irradiation 1000 W/m? and at the temperature values 25 °C, 50 °C and 75 °C. We also display, 
in Figure 6(b), the curves AE=f(V). It is clearly that the simulated values of the model are in harmony with 
the measured data. In the voltage variation interval [0, Vmp], the theoretical and experimental I=f(V) curves 
are in good agreement and the absolute current error tends to cancel. When, at a voltage variation in the 
interval [Vinp, Voc], the current decreases rapidly and the AE increases slightly but still remains below 0.2 for 
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the temperature 25 °C and 0.8 for the temperatures 50 °C and 75 °C. In the [Vmp,Voc] area the voltage is 
almost constant while the current is exponential and is sensitive to a few simple variations. 


9} G=1000w/m? | | 


o Experimental data 
1} |— Simulated data 4 


0 1 L 1 L 1 
0 5 10 15 20 25 30 35 


G=1000W/n? 


Figure 6. Accuracy of the proposed method (a) measured and simulated I(V) characteristics and (b) absolute 
error AE(V), for KC200GT panel under varying temperature 


4. CONCLUSION 

This article proposed an efficient hybrid method to determine the five PV module parameters using 
the measured current-voltage characteristics and the simple diode model. This method is based on the 
analytical calculation of the series resistance using three specific points of the experimental current-voltage 
curves. The saturation current, the photo current and the parallel resistance are expressed depending on the 
ideality factor by solving of a system of three linear equations using the linear least squares approach. The 
ideality factor value is determined iteratively to minimize the NRMSE value. The results obtained for two 
experimental cases studies of solar cell and PV panel show a good agreement and low error between the 
measured and the simulated data. 
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